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Abstract 

In this paper we formulate a geometric theory of thermal stresses. Given a temperature distribution, 
we associate a Riemannian material manifold to the body, with a metric that explicitly depends on the 
temperature distribution. A change of temperature corresponds to a change of the material metric. In 
this sense, a temperature change is a concrete example of the so-called referential evolutions. We also 
make a concrete connection between our geometric point of view and the multiplicative decomposition of 
deformation gradient into thermal and elastic parts. We study the stress- free temperature distributions 
of the finite-deformation theory using curvature tensor of the material manifold. We find the zero-stress 
temperature distributions in nonlinear elasticity. Given an equilibrium configuration, we show that a change 
of the material manifold, i.e. a change of the material metric will change the equilibrium configuration. In 
the case of a temperature change, this means that given an equilibrium configuration for a given temperature 
distribution, a change of temperature will change the equilibrium configuration. We obtain the explicit form 
of the governing partial differential equations for this equilibrium change. We also show that geometric 
linearization of the present nonlinear theory leads to governing equations that are identical to those of the 
classical linear theory of thermal stresses. 
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1 Introduction 

Classical elasticity theory quantifies the amount of stretch in a body by using a specific configuration as the 
reference configuration. The displacements as measured from the reference configuration and the strains asso- 
ciated with them are then used to get the stresses via constitutive relations. This viewpoint works nicely when 
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there is a relaxed, stress-free configuration that can be used as the reference configuration. However, this is not 
always the case. A body may have various sources of residual stresses, e.g. defects such as dislocations and 
disclinations, in which case there may not exist any stress-free configuration. One can observe the existence 
of these residual stresses by cutting pieces from the body when there are no external forces, and see that the 
pieces relax upon being cut. One can deal with these residual stresses in the classical theory [5] but we will 
use a different approach in this paper; that of geometric elasticity, and the notion of a space that describes the 
intrinsic properties of a body with a residual stress distribution. 

Temperature enters free energy density as a state variable. In classical linear theory of thermal stresses 
[5], it is assumed that there exists a reference temperature T at which the body is stress free. Free energy is 
then expanded about T and only linear and quadratic terms are kept. The governing equations of this theory 
consist of those of linearized elasticity and heat conduction with some coupling terms. Given an equilibrium 
configuration of the body at temperature T, a change in temperature will change the equilibrium configuration 
due to the coupling terms. Similar ideas are used in the nonlinear theory by looking at thermal stresses as a 
coupled nonlinear elasticity/heat conduction problem. 

In this paper, we study thermal stresses geometrically by considering a material manifold that explicitly 
depends on temperature. Material manifold is endowed with a Rlemanian metric G. We assume that given 
a reference temperature distribution To, when the body is unloaded in a Ricmannian manifold (B, Go) it 
is stress free. Change of temperature changes the metric. For similar ideas in the case of dislocations see 
[4; 15; 16; 17; 18; 20; 19]. We should emphasize that here we assume that temperature distribution is given. A 
geometric formulation of the coupled elasticity/heat conduction will be discussed in a future communication. 

As a motivation for our viewpoint, consider a piece from a thin, elastic spherical shell being forced to lie 
on a plane, e.g., by being squeezed between two flat surfaces. This constraint will induce stresses on the shell. 
Let us for the moment imagine that we are observing this shell from the two-dimensional viewpoint of the 
plane, ignoring the third dimension. When we cut pieces from the shell (all still forced to lie on the plane), we 
will observe that the former relax by a certain amount, demonstrating the "residual stress" on the body. The 
two-dimensional, planar viewpoint dictates that there is no stress-free configuration for this piece of material. 
However, if this same shell is "forced" to live on the surface of a sphere of appropriate radius, there will be a 
stress-free configuration, which may then be used as the reference configuration for measuring the amount of 
stretch, etc. The surface of the sphere is intrinsically two-dimensional, in the sense that it can be described 
without any reference to a third dimension, by using only two coordinates, and an intrinsic measure of distance, 
e.g. by using the spherical coordinates: 

ds 2 = R 2 (d0 2 + sin 2 6 d(/> 2 ). (1.1) 

This suggests a generalization to three dimensions: given a body with residual stresses with no apparent stress- 
free configuration, could it be possible to find some abstract three-dimensional space in which it would be 
stress-free? As stated, this is a very general question, and the answer depends on the source of the residual 
stresses for the case at hand and the notion of "space" one is considering. 

In this paper, we answer this question affirmatively for a specific source for residual stresses: that induced 
by the thermal expansion due to a non-uniform temperature distribution in a solid material. Our analysis is 
constructive; given a temperature distribution on a previously stress-free material, we construct a Riemannian 
manifold (to be defined shortly) in which the material under consideration with the given thermal profile 
would have a stress-free state. This construction is not of purely academic interest, since given the appropriate 
constitutive relations, it allows us to calculate the stresses in an elastic body in a given thermal setting, and is not 
restricted to a linear response. The appropriate way to quantify deformations for a body with residual stresses is 
by considering a map ( "the configuration map" ) from the "material manifold" to the three-dimensional spatial 
space that the body lives in. The constitutive relations are then written in terms of this map. 

In the remaining of the paper we study the effect of changes of the material manifold on equilibrium config- 
uration. We do this in a general setting when the material manifold is Riemannian. We also consider the case 
where the change in metric is small enough so that a linear approximation would be enough. In the case of small 
changes in material manifold we obtain the governing equations of evolution of the equilibrium configuration as 
a function of changes of the material manifold. 

Let us pause here and explain the main ideas and conclusions. Suppose we have a flat, two-dimensional 
clastic material on a plane. If we heat this material in a non-uniform way, it will tend to bend out from the 
plane, and take the shape of a curved surface. If we then force this curved surface to live in the flat plane by 
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perhaps squeezing it between two flat, rigid surfaces, we will induce some stresses on it. This suggests that for 
a given non-uniform temperature distribution, there may be a curved, stress-free shape that an "originally flat" 
elastic material wants to take. When we force this material to live on a flat plane, we induce "thermal stresses" . 
Of course in real life, we have three-dimensional elastic bodies, forever bound to live in flat, three-dimensional 
Euclidean space. 

Here is the main idea: given a material metric G describing the stress-free state of an elastic material 
at constant temperature, and a spatial metric g (for simplicity, both of these metrics can be taken as flat, 
three-dimensional metrics), we claim that a non- uniform temperature distribution on the material should be 
represented as a change in the material metric G. In particular, for a material with isotropic thermal properties, 
we claim that the required change in G is just a pointwise rescaling (i.e. multiplying the metric with a scalar 
function on the material manifold), and the way the scaling factor depends on the temperature is determined 
by the physical properties of the material. Non-isotropic thermal expansion will also be considered. 

We believe this is an example where a change in the material manifold can be clearly understood conceptually. 
This may be very helpful in understanding the role of the change in material connection in defect mechanics, 
where one has to consider not only a change in the metric, but also the "torsion part" of the connection. In 
particular, we have the hope that an analogy with thermal expansion will help clarify the often encountered (and, 
in our opinion, confusing) discussions of decomposing the deformation gradient F into "elastic" and "plastic" 
parts. The analogue of a plastic deformation in our case is a change in the temperature, resulting in a change 
in the material connection. Such a change in the material connection may induce stresses on a previously 
stress-free configuration, and as a result, change the equilibrium configuration. The resultant stresses in the 
new equilibrium configuration should be explored with the "usual" geometric elasticity. 

This paper is structured as follows. In §2 we motivate the connection between thermal stresses and changes 
in material metric. We do this by looking at the example of a two-dimensional disk and study the possibility of 
existence of a relaxed configuration embedded in the three-dimensional Euclidean space for an arbitrary radial 
temperature distribution. We also study the stress-free temperature distributions in nonlinear elasticity. In §3 
we present the main ideas of a geometric formulation of thermal stresses. We study the effect of a change of 
material manifold on the equilibrium configuration. In §4 the geometric linearization of the nonlinear theory is 
presented. Conclusions are given in §5. To make the paper self contained, we briefly review the basic concepts 
of diffcrcntoial geometry, parralclizable manifolds, and geometric theory of elasticity in the appendix. 

2 The Material Metric and Non-Uniform Temperature Distribu- 
tions 

Motivation. Suppose we start with a stress-free isotropic material with a uniform temperature distribution 
T\ and free boundary conditions, and increase its temperature to T2. The material will expand, and the original 
distance 5L\ between two neighboring points A and B in the solid body will increase to 8L2- The quantity 
(5L2 — SLi)/SLi turns out to be independent of the two points, i.e., the expansion is uniform. Note that 



where a is the coefficient of thermal expansion. Let us now assume that we use a Lagrangian coordinate 
system, X 1 , X 2 , X 3 , 1 !^., assume that the same material points have the same coordinates before and after the 
expansion. Then, the distance between the two points A and B is given approximately in terms of the metric 
tensor Ga as follows: 



where the components Gij are evaluated at a point between Xa and Xb- This shows that Gij should somehow 
depend on temperature. In other words, this suggests that, in this Lagrangian setting, we should be using 
different metric tensors for T\ and T2. Note that this relation between the Lagrangian coordinates and the 
material manifold works only because the material is in a relaxed, stress-free state in both temperatures. 
Otherwise, the distance that the material metric G measures would not be the spatial distance. 



SL2 — SLi 



a{T 2 - Ti) 



(2.1) 




(2.2) 



1 The superscripts denote coordinate labels. 
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The thermal expansion coefficient. Let us now connect the above description of the material metric 
in terms of a temperature distribution to the thermal expansion coefficient used in the classical theory. We 
imagine a material with a nonuniform temperature distribution, for various values of the constant temperature. 
The thermal expansion coefficient is defined by looking at the equilibrium volume of the material at different 
temperatures. Let us therefore look at the volume element of a material at a given constant temperature. We 
assume that the material metric is given by 

Gu(X,T) = H u (X)e 2 ^ T \ (2.3) 

where Hu is independent of temperature and T = T(X). The volume form associated with this metric is given 

as 

dV(X, T) = y/dct \Hu\ e N ^ T) d N X 1 (2.4) 

where d N X is shorthand for dX 1 A dX 2 • • • A dX N 2 Differentiating with respect to T, we obtain 

^dV(X,T) = VdetW e N «^d»X ™£1 = dV(X,T)N^l. (2.5) 

Thus, we can read-off the thermal expansion coefficient in terms of the temperature dependence of w(T): 3 

. , du>(T) 

a{T) =ar- (2 - 6) 



Remark: Suppose Hjj = 5jj 1 i.e. the initial material manifold is Euclidean. In this case G/j is conformally 
flat. It is known that any 2-dimensional Ricmannian manifold is conformally flat and the function lo is unique 
[2]. See the appendix for more discussions on this. 



2.1 Stress-Free Temperature Distributions 

Before we go into the general geometric theory of thermal stresses, we will first demonstrate an important 
application of our geometric approach, that of stress-free temperature distributions. 



The two-dimensional case. As in the introduction, consider a two-dimensional shell restricted to live on a 
flat planar surface between two rigid planes. We will assume that initially the shell is at constant temperature, 
and is stress-free, with no external or body forces. We would like to find the temperature distributions that will 
result in equilibrium configurations with zero stress. Changing the temperature uniformly will result in uniform 
expansion, and hence no stress. Are there other temperature distributions with this property? The answer is 
yes as is already known in the framework of linear thcrmoclasticity [5]. Due to the nature of our geometric 
approach, we will also be able to answer this question in a nonlinear setting. To the best of our knowledge, this 
has not been done in the literature. 

The spatial distances between material points are measured by the ambient space metric (the "spatial 
metric" ) , which is Euclidean. A given temperature distribution will result in a change in the material metric, as 
described above. A configuration will be stress-free if there is no "stretch" in the material, i.e., if the material 
distance between two points is the same as the spatial distance. This happens only if the two types of metric 
tensors (spatial and material) agree on the distance measurements between nearby material points, i.e. only if 
they are isometric. Since the spatial metric is assumed to be Euclidean, this means that the material metric, 
after the change due to a given thermal distribution, has to be Euclidean. 

It is worth emphasizing that one cannot simply set Gu = Sjj, the precise requirement is that the pull-back 
of the spatial (Euclidean) metric by the deformation map ip has to be equal to the material metric. This issue 
is closely related to the fact that a metric may be Euclidean "in disguise", i.e., one can write the flat two- 
dimensional metric in different coordinate systems, and it is not always easy to recognize that the metric is flat 
by simply looking at its components in a given coordinate system. 



2 Notc that N = 2 or 3. 

3 Note that for thermally isotropic materials volumetric thermal strain is N times the thermal strain. 
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Riemann's original work solves this problem for any dimensionality by defining the curvature tensor of the 
metric: a metric is flat, i.e., it can be brought into the Euclidean form 5u locally by a coordinate transformation, 
if and only if its curvature tensor is zero [2]. 4 

It turns out that in two dimensions, a weaker requirement is sufficient [2]: a metric is flat if and only if its 
scalar curvature (the Ricci scalar) is zero. Let us now apply this condition to a two-dimensional metric that is 
obtained from a non-uniform temperature distribution on an initially stress- free, planar shell, i.e., Gjj = e 2n Sjj, 
where O(X) = w(T(X)). The Ricci scalar for a metric of this form is given by [40] 

R = -2er 2n V 2 fl. (2.7) 

Thus, R = requires V 2 = 0, i.e., the exponent in the scale factor has to be a harmonic function. If we assume 
that ui(T) depends on temperature linearly, we obtain V 2 T = 0. This is exactly the same condition encountered 
in linear elasticity, see, e.g. Boley and Weiner [5]. This means that for the case of constant thermal expansion 
coefficient in two dimensions, harmonic temperature distributions do not result in any stresses. However, 
our result is more general: even if the thermal expansion is non-linear, we obtain the condition X7 2 fl(X) = 
\7 2 uj(T(X)) = 0, where u>(T) gives the general, non-linear dependence of (isotropic) thermal expansion to 
temperature. 

It is worth emphasizing the distinction between local and global flatness, and the implications for stress-free 
thermal distributions. Even though the surface of a right circular cylinder in three dimensions looks curved, 
it is locally, intrinsically flat. For any given point on the cylinder, one can find a finite-sized region containing 
the point, and a single- valued coordinate system on this region, for which the metric has the Euclidean form. 
Physically, this means that for any given point, we can cut some finite-sized piece containing the point, and can 
lay the piece on a flat plane, without stretching it. The surface of a sphere in three dimensions, on the other 
hand, is intrinsically curved; it is impossible to make any finite-sized piece of the sphere, no matter how small, 
to lie on a flat plane without stretching it. 

Curvature conditions like R = 0, or V 2 = can only detect such local issues. That it is impossible to make 
a full cylinder lie in a plane nicely (i.e., without tearing, folding, or stretching it) is due to the global topology of 
the cylinder, and local restrictions on curvature are not capable of constraining the global properties sufficiently 

In the context of thermal stresses, this subtlety is nicely demonstrated by the following example. Let us 
specialize to the case where Q depends only on the radial coordinate R of an initially flat annular piece of a 
material, Rq < R < R\. The flatness condition gives 

Solving this gives 

e 2 ° = 1 R 2P , (2.9) 

where 7 > and /3 are constants. 5 Thus, we are concerned with temperature distributions that result in metrics 
of the form 

dS 2 = R 2f} (dR 2 + R 2 dO 2 ) , (2.11) 

where we have set 7 = 1 by a rescaling of R. It may not be immediately obvious that these metrics are flat, 

1 

but a transformation to a new radial coordinate rbyi? = r' ) + 1 gives 

dS 2 =c 2 dr 2 +r 2 d0 2 , (2.12) 
where c = -g^j- A further transformation, R = \c\r, = 0/|c| gives 

dS 2 = dR 2 + R 2 d0 2 , (2.13) 



4 Wc have collected the definitions of various curvature-related quantities of interest in the appendix. 

5 Assuming T(Rq) = To, and that the coefficient of thermal expansion is a constant a, this solution corresponds to the temperature 
distribution 

T(fl)=T + ^ln(A) . (2.10) 
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which is the flat two-dimensional metric, except for an important subtlety. In the original coordinates, a point 
(R, 0) was identified with the point (R, + 2ir). In terms of the new coordinates, this means that a point 
(R, 0) needs to be identified with the point (i?, + 2tt/\c\). The geometric meaning of this is clear: the metric 
(2.13), with the proper identifications, is describing an annular piece from a conical surface, with deficit angle 
a = 2tt(1 - l/|c|), see Fig. 2.1. 



Figure 2.1: Zero-stress deformation of an annulus to a cone. 

Now, intuitively, one can guess that it will be impossible to make such a conical surface lie on the plane 
without tearing, stretching, or folding it. Thus, if we start with an annular shell between two rigid planes, a 
temperature distribution of the form (2.10) will indeed result in stresses, although the related material metric 
is intrinsically flat. However, if the material consists only of a simply-connected piece of the annulus (say, 
Ri < R < i?2, < 0i < < 02 < 27r), the temperature distribution (2.10) will just cause a stress-free 
expansion of the material, between the two rigid planes. See Fig. 2.2. 



Figure 2.2: Zero-stress deformation of a simply- connected piece of an annulus. 

This issue exists in the classical theory of elasticity, as well, see p. 256 in [5] where a set of global conditions 
have been given for stress-free thermal distributions in two-dimensions for multiply-connected bodies. As an 
application of these ideas, suppose we are given a temperature distribution T(X), and are seeking a temperature- 
dependent coefficient of thermal expansion a(T) that would result in a stress-free equilibrium state. In general, 
this amounts to solving the equation 

da dT dT „ dR d 2 T .„ da, ,o , 

This equation may or may not admit a solution, depending on T(X), but for the simple case of a radial 
temperature distribution T(R), the solution (2.9) dictates that a(T) is given by 

^ = % = % = Tmwy (2 - 15) 

an v ' 

As above, however, this works only for simply-connected pieces such as those described by 0! < < 2 and 
R\ < R < i?2, with 02 — 0i < 27r. With this in mind, a(T) for some specific radial temperature distributions 
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are as follows. 



T(R) 



T(R) = 
T X R X 



T Ri — TiR 



R\ — Ro 

TqRq 



Ti - To 



Ri — Rq 

In (r?/R^ 



Ri — Ro 
RqRi{To — Ti) 



]n(Rx/R ) 



R\ — Rq 

Tx-Tq 
In {Ri/R ) 



R 
1 

R 



a(R) = ctQ 
a(R) = Qo 



Ro 
~R' 
R 
Ro 



In R : a(R) = ag. 



(2.16) 
(2.17) 

(2.18) 



Let us next consider the case where the two-dimensional material is allowed to bend into the third dimension, 
instead of being squeezed between two rigid flat planes. Assume, once again, that we start with a stress-free, 
planar piece of material at uniform temperature, and we introduce a non-uniform temperature distribution. 
Will the material be able to find a stress-free state by bending into the third dimension? We will investigate 
this problem for the special case of a radial temperature distribution, and cylindrically symmetric configurations 
in three dimensions. 

According to our approach, the metric 



dS 2 



2a>(T(R)) 



(dX 2 + dY 2 



a 2u{T(R)) 



(dR 2 + R 2 d0 2 ), 



(2.19) 



describes the natural, stress-free state of the shell, and for a given configuration in the ambient 3-dimcnsional 
space, stresses will be due to a different metric being induced on the shell by the embedding in the ambient 
space. We will seek embeddings for which the metric induced from the ambient space is the same as the intrinsic 
metric, which, as opposed to the case considered above, is not necessarily flat. 

Let us begin by writing the induced metric for a given configuration with cylindrical symmetry, which is 
best done in cylindrical coordinates (p, (f>, z). Instead of the most general configuration, we will seek a solution 
of the form 



(j>(R,e) 
p(R,Q) 
z(R,G) 



e, 

p(R), 

z(R). 



For such an embedding, the metric induced from the ambient space is given as 

2 / 7 \ 2 



dS; 



2 

induced 



dR A 



dz 
dR 



dp 
dR 



p{R? 



(2.20) 
(2.21) 
(2.22) 



(2.23) 



In order for this induced metric to be the same as the intrinsic, material metric given by (2.19), we need 

p(R) = Re n{R \ (2.24) 



dz 
dR 



dp 
dR 



where we define fl(R) = tu(T(R)). Substituting p(R) in the second equation, we obtain 



^y + e ^(R) [l + i? ^ (i?) ] 2 = e 



2Cl(R) 



(2.25) 



(2.26) 



which, in principle, lets us solve for z(R). For a region where p(R) is invcrtiblc, we can also obtain the surface 
in three dimensions as given by z(p), by solving 



[i + Rn'(R( P ))Y 



(2.27) 



Note that for a uniform temperature distribution T(R) = To, £l'(R(p)) = and hence z{p) = zq, which is 
what we expect, i.e., in this case the relaxed configuration is planar. Note also that (2.27) has a solution only 
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if —2/R < Q'(i£) < 0. In order to have some insight about the meaning of this constraint, let us specialize 
to conical metrics described above in (2.13), e 2a = "fR 2f3 . For this case, the constraint —2/R< 0'(i?) < 
translates to —2 < j3 < 0, which gives, in terms of the deficit angle a, < a < 2ir. The condition a < 2ir is not 
surprising, but there is nothing wrong with a cone with negative deficit angle in terms of intrinsic geometry. 
The lower bound on f2 is simply telling us that it is not possible to embed such a cone in R 3 in a cylindrically 
symmetric way, which makes intuitive sense, considering the twisted shape of a saddle. 



The three-dimensional case. Let us next consider the three-dimensional case. In three dimensions, a 
vanishing Ricci scalar is not sufficient to guarantee local flatness, however, a three-dimensional metric is flat if 

„2fl/^(0) 



and only if its Ricci tensor vanishes [2]. The Ricci tensor TZjj of the metric Gu = e Gj j is given in terms of 
the Ricci tensor Tcfj of Gj] by the following relation [40] : 

IZu = Tlf) - (n- 2)V/VjO - Gf] (g (0) ) KL V k V l O + (n- 2)V/OV./0 - (n- 2)Gf) (g (0) ) AL V a -OV l O , 

(2.28) 

where n is the dimensionality. Now, once again, assume that the initial metric G*f} = 5u, — 0, and n = 3, 
and replace the covariant derivatives with partial derivatives. This gives 

Ku = -djdjQ - S IJ S KL d K d L n + d&djtl - S IJ S KL d K nd L n = 0. (2.29) 

Similar to what is implicitly done in classical linear thcrmoelasticity, let us assume that the reference tem- 
perature To is uniform, i.e. independent of position and that change of temperature is "small" . This means 
that Tjyr is small. But note that 

dn , N dT , 

aF = a ^ (2 ' 30) 

Therefore, <9/f2 is small too, i.e. quadratic terms in <9/f2 can be ignored. This gives us the condition that all 
the second derivatives of Q have to vanish. This means that Q is a linear function of the original Euclidean 
coordinates. If we further assume that w(T) is a linear function of temperature, we see that temperature itself 
has to be a linear function of the original Euclidean coordinates. Therefore, we recover the classical result 
in linearized thermal elasticity that in three dimensions, the only stress-free temperature distributions for an 
initially stress-free material depend linearly on the coordinates [5] . 

In the nonlinear case IZij — is equivalent to the following system of nonlinear partial differential equations 
in terms of Q: 

(2.31) 
(2.32) 
(2.33) 
(2.34) 
(2.35) 
(2.36) 

Let us first look at the following nonlinear partial differential equation for w = w(x,y). 

d 2 w dw dw 



0,12 


= 4 fi, 2 


) 






0,13 


= fi,lfi,3 


) 






0,23 


= 0, 2 0, 3 


) 






0,n 


+ V 2 Q + 


2 2 


+ 2 3 


= o, 


0,22 


+ V 2 Q + 


Ofi 


+ 2 3 


= 0, 


0,33 


+ V 2 Q + 


2 i 


+ 2 2 


= 0. 



dxdy dx dy 
Using the change of variable u = e~ w , one obtains 

d 2 u 



(2.37) 



dxdy 

Thus, the most general solution is 



= 0. (2.38) 



(x,y) = -Hf(x)+g(y)], (2.39) 



2.2 Connection with Multiplicative Decomposition of Deformation Gradient 9 



for some arbitrary functions / and g. Using (2.31)-(2.33), one can show that 

Q{X\ X 2 , X 3 ) = - HfiX 1 ) + g{X 2 ) + h(X 3 )], 
for some arbitrary functions f,g, and h. Now substituting (2.40) into (2.33)-(2.36), we obtain 

f"(X 1 )=g"(X 2 ) = h"(X 3 ). 

Therefore 

f(X 1 )=c (X 1 ) 2 + d 1 X 1 + d 2 , 
g(X 2 ) = c (X 2 ) 2 + d 3 X 2 + d 4 , 
h(X 3 ) = c (X 3 ) 2 + d 5 X 3 + d 6 , 

for some constants Co, di, c?6. 6 Plugging these back into (2.34)-(2.36), we get 

3 



In 



i=l 



Shifting the origin X % — > X 1 + b 1 , this becomes 

n(x\x 2 ,x 3 ) = - In (c R 2 ) , 

where R = ^{X 1 ) 2 + (X 2 ) 2 + (X 3 ) 2 . 

If a is constant, then this corresponds to the following temperature distribution 7 



(2.40) 



(2.41) 



(2.42) 



(2.43) 



(2.44) 



In (c R 2 



a{T - T ) or T - T = c 



lni?. 



(2.45) 



Note that Co represents a uniform change in temperature. In order to understand what this solution represents 
physically let us write the metric in polar coordinates. 



dS 2 = e 
Now let us define 



2! 2 



1 



[dR 2 + R 2 (de 2 + sin 2 Qd® 2 )] = [dR 2 + R 2 (d0 2 + sin 2 9d$ 2 )] 



In terms of R, the metric becomes 



R= — . 

cR 



dS 2 = dR 2 + R 2 (d0 2 + sin 2 6d</> 2 ), 



(2.46) 
(2.47) 

(2.48) 



which is precisely the flat Euclidean metric in three dimensions. Thus, after the thermal expansion, the metric 
is still flat, but the radial coordinate in which it is manifestly so is related to the old radial coordinate by (2.47) 
(up to a simple shift of origin). This means that, particles at the two radii R\ < i?2 move to the new radii 
R\ > i?2, after the thermal expansion, i.e., the material gets "inverted". This may not be possible for a solid 
ball without tearing it apart, but it is perfectly possible for a piece from such a ball, as demonstrated in Fig. 2.3. 



2.2 Connection with Multiplicative Decomposition of Deformation Gradient 

Geometric study of thermal stresses goes back to the works of Stojanovic, et al. [37]; Stojanovic [38]. These 
researchers extended Kondo's [15; 16; 17; 18] and Bilby's [4; 3] idea of local clastic relaxation in the continuum 
theory of distributed defects to the case of thermal stresses. See also Maugin [32] for a review of some relevant 

6 Note that the case cq = di = d% = d§ = corresponds to uniform temperature distributions. 
7 Note that this is similar in form to the 2D solution (2.18). 
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Figure 2.3: A nonlinear stress-free thermal deformation of a ball. Note that because of symmetry only deformation of a great half 
circle of the projected sphere on a plane passing thorugh the center of the ball is shown in this figure. 



works. One should note that the idea of using differential geometry in anelasticity goes back to an earlier work 
by Eckart [7]. 

Stojanovic's idea is similar in spirit to the approach described in this paper: a nonuniform temperature dis- 
tribution in general leads to residual stresses essentially because the body is constrained to deform in Euclidean 
space. If one partitions the body into small pieces, each piece will individually relax, but it is impossible to re- 
alize a relaxed state for the full body by combining these pieces in Euclidean space. Any attempt to reconstruct 
the full body by sticking the particles together will induce deformations on them, and will result in stresses. An 
imaginary relaxed configuration for the full body is incompatible with the geometry of Euclidean space. 

The approach taken in this paper is to ask the question: which space, as opposed to the Euclidean space, 
would be compatible with a relaxed state of the body? We claimed above that the answer to this question is, a 
Ricmannian manifold whose metric is related to the nonuniform temperature distribution by (2.3). This metric 
describes the relaxed state of the material, and the strains in a given configuration should be measured with 
respect to the relaxed state, i.e., the new material metric. In this setup, the constitutive relation (e.g., a free 
energy function) is given in terms of the material metric, the (Euclidean) spatial metric, and the deformation 
gradient F. The constitutive relation allows one to calculate the stresses induced for a given configuration, at 
least in principle. 

Stojanovic, on the other hand, takes the following viewpoint. Consider one of the imaginary relaxed pieces 
described above. The process of relaxation after the piece is cut corresponds to a linear deformation of this piece 
(linear, since the piece is small) 8 . Let us call this deformation Fy. If this piece is deformed in some arbitrary 
way after the relaxation, one can calculate the induced stresses by using the tangent map of this deformation 
mapping in the constitutive relation. 

Now, in order to calculate the stresses induced for a given deformation of the full body, we focus our attention 
to one such particular piece. The deformation gradient of the full body at this piece F can be decomposed as 
F = F e Fy, where, by definition, F e = FF^ 1 . Thus, as far as this piece is concerned, the deformation of the full 
body consists of a relaxation, followed by a linear deformation given by F e . The stresses induced on this piece, 
for an arbitrary deformation of the full body, can be calculated by substituting F e in the constitutive relation. 

Note that F e and Ft are not necessarily true deformation gradients in the sense that one cannot necessarily 
find deformations ip e and tpx whose tangent maps are given by F e and Ft, respectively. This is due precisely 
to the incompatibility mentioned above. However, as long as we have a prescription for obtaining F e and Ft 
directly for a given deformation map if for the body and a temperature distribution, we can calculate the stresses 
by the following prescription. 

For an isotropic material, Stojanovic gives the following formula for Ft in terms of the temperature. 

(F T ) A B = §{T)5g. (2.49) 

This means that a small piece relaxes by a uniform expansion, whose magnitude is determined by a function $(T) 
that characterizes the thermal expansion properties of the material under consideration. Given this formula 
for Ft, we can calculate F e = FF^ 1 for a given deformation, and utilize a constitutive relation that gives the 
stresses in terms of F e . 

These two approaches seem very different philosophically, and at first sight, Ft, the "incompatible interme- 
diate deformation gradient," perhaps seems a little mysterious from the geometric standpoint. However, these 
approaches are related, as we will demonstrate next. Our discussion can easily be generalized to other sources 
of residual stresses; a local relaxation approach and the Ricmannian approach are equivalent for a large class of 



This transformation is not necessarily uniquely determined, see below for a discussion of this issue. 
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settings. 9 

As mentioned above, Stojanovic [37; 38] gives (Ft) A b = "&{T)5g and relates the coefficient of thermal 
expansion a to $(T) by 

oc{T)=d{T)^-. (2.50) 



This can agree with (2.6) only if 



e 



,(T) 



•&{T). (2.51) 

In order to show the mechanical equivalence of the two approaches by using this identification, we need to show 
that for any given constitutive relation for one of the approaches, one can find a corresponding constitutive 
relation for the other approach that predicts the same stresses for all possible deformations when § and ui are 
related through (2.51). 

The constitutive relations of the two approaches are formulated in terms of different quantities: F e = FF^ 1 
on one side, and G(T) and F on the other. Let us start with our approach, namely, assume that a constitutive 
relation is given in terms of G(T) and F. This takes the form of a scalar free energy function that depends on 
G(T), F, as well as on the spatial metric tensor g, and possibly X and T(X) explicitly: 

* = *(X,T,G(X,T),F,g). (2.52) 

Now, G, F, and g are tensors, written in terms of specific bases for the material space and the ambient 
space. Commonly, bases associated to coordinate systems are used. A change of basis changes the components 
of these tensors, but 'J, being a scalar, does not change. Let us consider a change of basis from the original 
coordinate basis E^ of the material space, satisfying 



(E a ,E s ))g = G AB , (2.53) 



to an orthonormal basis E^ that satisfies 



({E A ,E B )) G =S A§ . (2.54) 
The transformation between the two bases is given by a matrix F A B as 

E A = F A B E B . (2.55) 

The orthonormality condition gives 

F/F/G CD =^. (2.56) 
Any F A G that satisfies this equation gives an orthonormal basis. Given such an F A C ', we can also obtain an 
orthonormal basis for the dual space by using its inverse. Defining F c o as the inverse of the matrix F A B , i.e., 
F A B F A c = S B and F a b F C b = S A , we obtain the dual orthonormal basis {E A } in terms of the original dual 
basis {E A } by 

E A = F i s E B . (2.57) 

For thermal stresses, assuming that the initial material manifold is Euclidean, Geo = e 2u ^5cD — ${T) 2 5cd 
gives 

F A C = 8 C A e-" (T) = 5 C A iT^T), (2.58) 

as a solution to (2.56). Here, S B is 1 for A = B, and 0, otherwise, i.e., 5^ = 5? = 6? = 1, etc. Note that (2.56) 
has other solutions, too, which we will comment on below. 

Now let us write the components of the total deformation gradient F in the orthonormal basis {E^}. The 
components are transformed by using F: 

pa A = ^A Bpa B ■ (2.59) 



9 However, there are cases that will require a further generalization, namely, cases where the material manifold has a connection 
with torsion and non-metricity. 
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Now, using (2.58), (2.51), and (2.49), we see that the components F a A are given precisely by those of F e , the 
"elastic part" of the deformation gradient in Stojanovitch's approach: 

pa A = FA Bpa B = Sfe'^F a B = ($(T))~ 1 5 B F a B = {F^) A B F a B = (F e ) a A . (2.60) 

Thus, Stojanovic's F e is nothing but the original deformation gradient, written in terms of an orthonormal 
basis in the material space. In passing, we have also shown that there is no need for a mysterious "intermediate 
configuration" as the target space of Ft, the latter just gives an orthonormal frame in the material manifold, 
and as such, can be treated as a linear map from the tangent space of the material manifold to itself. As 
mentioned above, these ideas can be generalized to other problems with residual stresses. 

Rewriting the constitutive relation (2.52) by using an orthonormal basis for the material manifold, we obtain 

* = *(X,T,G AB = 5 AB ,F a B = (F e ) a B ,g ab ). (2.61) 

Thus, given a constitutive relation \I/ Rlem in the Riemannian approach, one can obtain a constitutive relation 
^ LR in the "local relaxation" approach by simply going to an orthonormal basis by (2.55) and (2.56), and 
ignoring the constant terms G AB = S A b and g a b = 5 a b in the functional dependence. 

* LR (X, T, (F e ) a B ) = * Ricm (X, T, G AB = S AB , F a B = (F e ) a B , 9ab = S ab ) . (2.62) 

Going in the opposite direction is also possible; starting with a free energy function for the Stojanovic's approach, 
one can derive an equivalent free energy in the Riemannian approach. This direction may be slightly more 
confusing, since the metrics of the material manifold and the spatial manifold are not explicitly written out 
initially. One proceeds by first writing F e in terms of its proper index structure (F e ) a b in the Riemannian 
approach, and inserting 5 a b and S A b where necessary for tensorial consistency, and finally interpreting these as 
the components of metric tensors, and performing a change of basis, if desired. 

Non-coordinate bases and torsion. Although a coordinate basis {E^ = d/dX A } is not necessarily or- 
thonormal, one can always obtain an orthonormal basis by applying a pointwise change of basis F A B . Moreover, 
giving an orthonormal basis in this way is equivalent to giving a metric tensor at each point; the inner product 
of any two vectors can be calculated by using their components in the orthonormal basis. We have seen above 
that in the context of thermo-elasticity, this means that a change in the material metric due to a change in 
temperature can be given in terms of the "thermal deformation gradient" of the local relaxation approach. 
Given an orthonormal basis {E^}, it is possible to obtain another one, {E^}, by using an orthogonal 

transformation A. A B as 

E^ = A/ E b , (2.63) 

where A A B satisfies A A C A^S^p = 5 AB - Let the relation between the original coordinate basis {F A } and the 
new orthonormal basis be given by the matrix f' A B as follows 

E^ = F' A B E B . (2.64) 

F' A B =A A d F d B . (2.65) 

Going in the opposite direction, one can see that F and F' represent the same material metric G, if and only if 
they are related through (2.65) for some orthogonal matrix A A B . This means that there is an SO(3) ambiguity 
in the choice of F, and hence, in that of Ft- 

As opposed to a coordinate basis, the elements of an orthonormal basis do not necessarily commute with 
each other; whereas [E A , E#] = for E A = d/dX A , for an orhonormal basis E^ = F a b Eb, one has [35] 

(2.66) 

where 

r . C _ pC (r E ^F B p E ^A \ (0 R7\ 



The relation between F and F' is given as 
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(J 

The connection coefficients V ^ B f° r an orthonormal basis, defined through 



= T% % (2.68) 



are related to the connection coefficients Y AB oi the coordinate basis by 

-C r DrC f dF i 



In a coordinate basis, the components of the torsion tensor are given by the antisymmctrization of the two 
lower indices of the connection coefficients, i.e. (see the appendix) 

Tab C = T C AB - T C BA . (2.70) 

However, for a non-coordinate basis, the components are given by 

ab = Tab ~~ ^ba ~ c ab^ ■ (2-71) 

Our formalism is based on Ricmannian geometry, and in particular, on the torsion-free Levi-Civita connection 
defined by the metric G(T). Thus, the torsion tensor for the material connection vanishes in both the coordinate 
basis, and the orthonormal basis. Let us show this explicitly for F A B = r d~ 1 {T)8 B ^. In the original coordinate 

basis, the metric tensor is given by Gab = s 2 ^ t ^Sab = $ 2 (T)Sab- Thus, the connection coefficients in this 
basis are 

Tsc = {#,bS£ + #, C S A - $ :D 5 AD 5 B c) ■ (2.72) 
Using (2.70), we have T A BC = 0. Next, using (2.67) and (2.69) with 

F A B = #6& F A B =^- 1 Sf, (2.73) 



we obtain 
and 



c AB 6 = ^ ~ - (2-74) 



^AB = -^-S%8l + l5ppSr DE C. (2.75) 

Using these in (2.71), we obtain T B( j A = 0. Note that the Riemann curvature tensor has the following form 

K A bcd = ?T 2 [2 {^ c Si - 0,u<$) 6,d + 2 {tf, B 6cD - V.cSbd) $,e8 ae + #,e$,f6 ef {S^Sbd - S a S C d)] 

+iT 1 [#,bdS$ - $,cdS a + {$,ce&bd - V.bdScd) S ae ] , (2.76) 

which, in general, does not vanish. 

Stojanovic, on the other hand, calculates a non-vanishing torsion tensor. This discrepancy is due to the fact 
that he uses the following connection 

r i (v-^\ A d(JY) M c (0 77 , 

r BC = ( F T ) M ~QXB ' ( 277 ) 

Similar connections have been used in other contexts [15; 4]. It can be shown that this connection has vanishing 
curvature but has nonvanishing torsion. This is related to the so-called canonical connection in absolutely 
parallclizable manifolds [8; 9; 44]. See also Epstein and Elzanowski [10] for similar connections in the context 
of inhomogeneities and their geometric representations. In the appendix, we give some details on absolutely 
parallelizable manifolds and the above connection. 

In summary, our geometric approach has a concrete connection with that of Stojanovic: in our approach 
we use a Ricmannian manifold with a temperature-dependent metric as the material manifold while Stojanovic 
implicitly uses the same metric but in an absolutely parallclizable manifold that is not Ricmannian. Using cither 
approach would be fine and a matter of taste, however we believe that our approach is more straightforward 
as we do not introduce an unnecessary torsion in the material manifold. Representing changes of temperature 
by a one-parameter family of conformal Ricmannian metrics enables us to find the zero-stress temperature 
distributions even in finite deformations. 
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2.3 Anisotropic Thermal Expansion 

So far, we have assumed that thermal expansion is isotropic, i.e. a change of temperature results in a change of 
length independent of orientation. Let us now see how one should modify the theory when a temperature change 
results in different changes in length in different directions and in possibly a change of shape. Even in the case 
of anisotropic thermal expansion all one needs is a temperature-dependent material metric G(X, T) but in this 
case the material metric is no longer a simple rescaling of the original material metric. The physical idea is the 
following. Given any point in the initial stress-free material manifold with metric Go(X), there exists a frame 
field {E 1 (X), E 2 (X), E 3 (X)} for each material point such that in this frame the material metric is diagonal and 
has the following form 

G(X, T) = e^^'E 1 <g> E 1 + e 2 " 2(T) E 2 ® E 2 + e 2 " 3 ( T )E 3 ® E 3 . (2.78) 
Given a coordinate basis Er = d/dX 1 , we have 

E/ = A J T Ej. (2.79) 

Thus 

G(X, T)=Y^ e 2u " {T) A J /Ej ® A K I 'E K . (2.80) 
i 

3 Geometric Elasticity with Temperature Changes 

Material manifold. As mentioned in the previous sections, the material manifold describes the intrinsic 
"shape" of the natural, stress-free state of the material. The geometry induced by a given configuration of the 
material in the spatial manifold may or may not agree with the intrinsic geometry The discrepancy between 
the two geometries (the induced and the intrinsic) is in general a cause for stresses, which are described by 
geometric constitutive relations. In this section, we will describe this framework. 10 

The motion of an elastic body is described by a possibly non-isometric, time dependent embedding of the 
material manifold in the spatial manifold (see Fig. 3.1). There is, however, another possible source of time- 
dependency: the geometry of the material manifold itself may change in time. 11 A change in the geometry of 
the material manifold is sometimes known as a referential change (see [12; 31] and references therein), and the 
precise meaning of this has sometimes been a source of confusion in the literature. In this paper, we have a case 
where the change in the geometry of the material manifold is described explicitly in terms of the temperature, 
and we believe that the conceptual clarity brought by this simple example may provide insights to other cases 
of referential changes, such as those that describe the evolution of defects in a crystalline solid. 




Figure 3.1: Motion of a continuum with temperature changes. 

10 Due to our approach to thermal stresses, in this paper, we treat the material and the spatial spaces as Riemannian manifolds, as 
in [29; 42], and by "geometry", we understand the Levi-Civita connection associated to the metric tensor. In general, the geometry 
of either of these spaces can be given by a more general connection that has torsion and / or non-metricity. Such connections have 
found use in the literature of defect mechanics. 

11 The geometry of the spatial manifold may also change, but we do not consider this issue in this paper. 
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In Section 2, we proposed to describe the thermal expansion of an isotropic material by a change in the 
material metric, given by equation (2.3): 

Gu(X,T) = H u (X)e 2 ^ T K (3.1) 

Here, lu(T) is a function that describes the thermal expansion properties of the isotropic material under consid- 
eration. The coefficient of thermal expansion is given by equation (2.6): 

aCT) = ^. (3-2) 

We assume that the temperature distribution in the material is given. If the temperature depends on time, then 
the material metric describing the relaxed state of the material will also depend on time, through (3.1). 

We should mention that evolution of reference configuration in the literature of continuum mechanics is 
more or less ambiguous. It is believed that an evolving reference configuration can model dynamics of defects. 
However, to our best knowledge, there are no concrete examples in the literature. We believe that the present 
geometric formulation of thermal stresses in terms of a temperature-dependent material manifold can make 
the role of reference manifold clearer and can shed light on other more complicated problems, e.g., continuum 
theory of solids with distributed dislocations. 

Conservation of mass. Let us begin by writing the conservation of mass in this setting. If the temperature is 
time- independent, the usual material version of the conservation of mass holds: the material density is constant. 

p (X,t) = p (X). (3.3) 

If, however, temperature changes in time, the material metric will expand or contract, so the material mass 
density will change. The evolution of the mass density will then be given by 

p (X,T)dV(X,T) = m (X), (3.4) 

where dV(X,T(t)) is the volume form of the metric G(X,T), and m(X) is the temperature-independent (and 
hence time- independent) differential form representing the mass density (mass form). This equation tells us 
that if the material manifold expands due to a temperature change, the total mass in a material region will not 
change, and hence the density po will decrease inversely with the increase in the volume of that region. Since 
the volume form is given by 

dV{X, T) = y/ dot \H U \ e N ^ T) d N X, (3.5) 

we can get the density for a given temperature T in terms of the density at an initial temperature To by using 
Po (X,T)dV(X,T) = p (X,T )dV(X,T ). This gives 

MX, T) = e N ^™-^ T » Po (X, T ). (3.6) 

In terms of the coefficient of thermal expansion a = 4^, this can be written as 



Po 



(X,T) = Po (X,To) e " JV/ -o Q ^. (3.7) 



Incompressibility. Elastic incomprcssibility means that elastic deformations cannot cause any changes in 
volume. Thus, for a given temperature distribution, the deformation map must preserve the volume element. 
The volume elements in the material and spatial manifolds, dV(X) and dv(x) are related by 

d«(x(X)) = J(X, T)dV(X, T), (3.8) 

where the Jacobian J is given as 

WT^detF^^. (3.9) 

Thus, incompressibility means that J(X,T) = 1. Given two Ricmannian manifolds, distance preserving maps, 
i.e., isometries between them may or may not exist. A similar question may arise for volume-preserving maps: 
given two Ricmannian manifolds (in our case, the material manifold and the spatial manifold), does there exist 
a volume-preserving map between them? Moscr [34] answers this question in the affirmative, so the study of 
incomprcssibility in this setting is not vacuous. 
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Free energy. The free energy, in addition to explicitly depending on temperature, will depend on the 
temperature-dependent material metric tensor as well, i.e. 

tt = *(X,T > G(X > T) J F,g). (3.10) 

Therefor, the first Piola-Kirchoff stress, given by 

P = P(X,T) = g- 1 — , (3.11) 
explicitly depends on the temperature-dependent material metric. 



Balance of Linear Momentum. Let us now look at the governing equations for a given temperature dis- 
tribution T = T(X). We will only study the static case 12 , for which the balance of linear momentum reads 

apaA 

DivP = Or p^ ]A = __ +V A B paB +1 a cF c A p b A = ^ (3^ 

where are the connection coefficients for the material metric Gab-, and are the connection coefficients 
for the metric g ab . This is the standard balance of momentum in geometric elastostatics, see, e.g., Yavari, et 
al. [42]. For the case of thermal stresses, the material connection coefficients T^ B are those of the metric (2.3), 
Gjj(X, T) = i//j(X)e 2w ' T ' ; they are given in terms of the connection coefficients ° °f the metric Hjj as 

[40] 

r % = r i A H B )c + (^8 B a + s^d A n - H AB H CD d D n) . (3.13) 

Suppose the initial material metric Hab is Euclidean. Then, using Cartesian coordinates, we have 

T% = 3d B n. (3.14) 
For a Euclidean spatial metric in Cartesian coordinates, g a b = 5 a b, we obtain 

In terms of the thermal expansion coefficient a = , this becomes 

QpaA dT 



3 a l^H P = °- ( 3 - 16 ) 



dX A dX 

In the following example, we show that in the geometric framework, some nonlinear problems can be solved 
analytically. 



Example. Let us consider a two-dimensional, incompressible neo-Hookcan material in a flat two-dimensional 
spatial manifold. The free energy density of a neo-Hookean material in two dimensions has the form 

* = *(X,C)=/i(trC-2), (3.17) 

where C is the Cauchy-Green tensor, or equivalcntly, the pull-back of the spatial metric, Cab = F a AF b Bg a b, 
and fx is a material constant. We will assume that this form holds for an isotropic material under thermal 
expansion, and in particular, we assume that there is no explicit temperature dependence in the free energy 
apart from the dependence through C. In components 

* = H (F a A F b B g ab G AB - 2) . (3.18) 

The "2" is of no particular significance: when the material metric is fixed, it just shifts the free energy by a 
constant. When the material metric changes as in (2.3), its contribution to the free energy is proportional to 



12 The governing equations for the dynamics case are similar. However, for the dynamic problem one has to consider an evolving 
temperature distribution governed by the heat equation. This will be discussed in a future communication. 
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the temperature-dependent material volume, which, for a given temperature distribution, is independent of the 
spatial configuration. We ignore this term, and use 'J = ^(X, C) = /itr C as our definition of the free energy. 

Let us assume that initially the material has a fiat annular shape R\ < R < R2 without any stresses, at a 
uniform temperature To. We would like to calculate the stresses that occur in the new equilibrium configuration 
after we change the temperature in a rotationally symmetric way, T = T(R). In polar coordinates, the spatial 
metric and its inverse read 

-(££)-(; s)-'t i -(££)-(; 1 aO- 

and thus dctg = r 2 . The only nonzero connection coefficients are: 

1'ee = ~r, lie = ll = 1/r. (3.20) 
For the temperature-dependent material metric we have 

G = ( ^ RR G J' B ) = ( I ^) e 2 -™), G 1 =( °!Z TZ \ = ( I . L ^ e ~ 2 "™\ 



G eR Gee J ~ \ R 2 J" ' ~ ~ V G eR G ee J~ \ 1/R 2 

(3.21) 

and thus, det G = R 2 e^ T( - R ^. The following nonzero connection coefficients arc needed in the balance of linear 
momentum: 

T R R = n'(R), T R e =-R- R 2 n'(R), r| e = r| fl = 1/R + n'(R). (3.22) 
In terms of the thermal expansion coefficient, these are given as 

T R R = aT'(R), rf e = -R- R 2 aT'(R), T% e = T% R = 1/R + aT'(R). (3.23) 

Given the temperature distribution T = T[R), we are looking for solutions of the form 

tp(R,Q) = (r,6) = (r(R),Q). (3.24) 

Thus 

T {)• F - - ( lf T ] ! ) < 3 - 25 > 



This gives the Jacobian as 

rr' 



(3.26) 



Re 2u(T) ■ 

Incompressibility dictates that 

rr' = Re 2 ^ T \ (3.27) 
This differential equation has the following solution 

r R 

r 2 (R)= 2ie 2 ^ T ^di + r 2 {R). (3.28) 
Jr.! 

Note that r% (R) is not known a priori and will be obtained after imposing the traction boundary conditions at 
r% and f2. 

In incompressible elasticity, P aA is replaced by P aA — Jp(F~ 1 )~ A i ) g ab , where p is an unknown scalar field 
(pressure) that will be determined using the constraint J — 1 [29], i.e. 

P aA = 2fiF a B G AB ^p(R)(F~ 1 ) A b g ab . (3.29) 

Therefore, using (3.27), we get the nonzero stress components as 

prR= 2 ! iR_ p{R) L e -MT(R)) > P ee = _ £^ (3 3Q) 

r R R z r 2 

where p(R) is an unknown pressure. 
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Balance of linear momentum in components reads 



P aA U 



For the radial direction, a = r, we have 



dP 



a A 



dX A 



Ti B P aB + P bA r b cF c A = 0. 



(3.31) 



f)prA 

P rA \ , — _l r A P rB a- P hA ^, r P c . 

P \A— K^Ta + 1 ABP + P lbc b A 



dX J 

QprR 

OR 

QprR 

dR 



n R +r% R )p rR ■ ^ 



P lee b e 



R 



2aT'(R) ) P rR - rP 



ee 



This gives 



p'(R) = 



2[iR 



2u(T(R)) 



R 



2 (1 + aRT') - ±^ e 2 ^™ - i- e -MT(«)) 



Assuming that p{Ri) = 0, we obtain 

-R 



p(R) 



2 M£ 2u(T(£)) 



2(1 + 0(0^(0) 



R? 



r 2 (0 



For a = 6, balance of momentum (3.31) gives 



96 



T A eP ee + p e Rl e rF r R + p-e^^ = (r « Q + r g e) p ,e = q 



(3.32) 



(3.33) 



(3.34) 



(3.35) 



i.e. this equilibrium equation is trivially satisfied. Therefore, given the temperature distribution T(R), we can 
calculate all the thermal stresses analytically. 



4 Linearized Theory of Thermal Stresses 

In this section, we linearize the governing equations of the nonlinear theory presented in the previous section 
about a reference motion. Geometric linearization of elasticity was first introduced by Marsden and Hughes [29] 
and was further developed by Yavari and Ozakin [43]. See also [33] for similar discussions. Here, we start with 
a temperature-dependent material manifold and its motion in an ambient space. Given a reference motion, we 
are interested in obtaining the linearized governing equations with respect to this motion. We will assume that 
the ambient space manifold is Euclidean. This is not a necessary assumption but it provides a natural setting 
for most practical problems of interest and will simplify the subsequent calculations. For simplicity, we will 
restrict attention to time-independent solutions. 

Suppose a given material with a temperature distribution T(X) and the related material metric G is in a 
static equilibrium configuration, ip. The balance of linear momentum for this material body reads 



Div P + p B = . 



(4.1) 



Now suppose we change the temperature of this material by a small amount <5T(X). This will change the 
material metric to G' = G + 5G, and ip will no longer describe a static equilibrium configuration. A nearby 
equilibrium configuration may be given by tp' = (p + 5ip, and the stress in this new equilibrium configuration 
will be P' = P + 6P. One would like to calculate the change in the stress (or the configuration), for a given 
small change in temperature (see Fig. 4.1). 

While the spirit of this setup is familiar from other linearization problems, some care is needed in interpreting 
its meaning. P(X) is a two-point tensor (it has components in both the material and the ambient spaces: P aA } 
based at X and ¥?(X), whereas P'(X) is based at X and p'(X). Defining <5P(X) = P' — P is nontrivial for a 
general ambient space metric. This is related to the fact that subtracting tangent vectors at different points 
in a manifold is only defined with respect to a choice of a path connecting the two points, on which a parallel 
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(B, G+5G) 



Figure 4.1: Motion of a continuum with temperature changes. 



transport is to be performed. By restricting our attention to a Euclidean ambient space we sidestep this issue, 
using the natural, path-independent parallel transport in Euclidean space. Another issue is the definition of 
§(p = ip' — p. While one can use coordinate systems to make approximate sense of this equation for two nearby 
maps, it is a little troublesome geometrically, since the subtraction of two maps between manifolds is not defined 
geometrically. 

The linearization procedure can be put to firmer footing if instead of talking about two nearby configurations 
and the differences of various quantities for these configurations, we describe the situation in terms of a 1- 
parameter family of configurations around a reference configuration, and calculate the derivatives of various 
quantities with respect to the parameter. These derivatives will capture the behavior of the solution as a 
function of the parameter, for small values of the latter. Thus, let T e (X) be a 1-parameter family of temperature 
distributions on our material manifold, G e be the corresponding family of material metrics, p e be the equilibrium 
configurations, and P e be the stresses. Let e = describe the reference equilibrium configuration. Now, for a 
fixed point X in the material manifold, <^ £ (X) describes a curve in the spatial manifold, and its derivative at 
e = gives a vector U(X) at ^(X): 

dp e (X.) 



U(X) 



(k 



e=0 



(4.2) 



Considering 8p> ~ e ^r~> we see that a more rigorous version of Sip is the vector field U. Similarly, one has 



SG 



de 



G f . 



e=0 



When the change in G is due to a change in T, we have 

du) dT 



d 
~dc 



e=0 



dT de 



6=0 



G = /3G, 



(4.3) 



(4.4) 



where /3 = 2^f| e=0 = 2a(T e )f| e=0 . 

Now consider, in the absence of body forces, the equilibrium equation Div P = for the family of temperature 
distributions parametrized by e: 

Div £ P e = 0. (4.5) 

Linearization of (4.5) is defined as [29; 43]: 



d 
d~c 



c=0 



(Div £ P e ) =0. 



(4.6) 



Once again, one should note that since the equilibrium configuration is different for each e, P e is based at 
different points in the ambient space for different values of e, and in order to calculate the derivative with 
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respect to e, one in general needs to use the connection (parallel transport) in the ambient space. For the 
Euclidean case we are considering and a Cartesian coordinate system x a , this is trivial. In components (4.6) 
reads 

dP aA (e) 



dX A 

Thus, the linearized balance of linear momentum reads 



T A B (e)P aB (e)=0 



d d 



dX A de 



=o 



P aA (e) 



e=0 



Note that 



paA _ gCic 



paB , -pA u 

F +lAB d~e 



<9* 



e=0 



P aB (e) = 0. 



(4.7) 



(4.8) 



dF c 



where $ = ^(X, T, F, G, g) is the material free energy density. In calculating 
changes in F and G due to the change in the equilibrium configuration: 



dP aA (e) 



dP aA (e) dP aA dF b B dP aA dG CD 



Defining 



A aA B 



QpaJ 



de dF h B de 

d 2 ^ 



dG 



CD 



de 



(4.9) 

we need to consider the 
(4.10) 



dF b B 9 dF b b9F c a " Geo 8Gcd9F c a 

where the derivatives are to be evaluated at the reference configuration e = 0. Noting that 



and B aACD 



jaA 



d 2 ^ 



dF a f 



de 



c=0 



8U a 
dX A '' 



wc obtain 



Using 



d_ 



P aA (e) = A aA b B U b 



e=0 



B 



aACD 



GcD (3- 



BC 



1 r<AD 
2° 



(dG 



BD 



dG 



CD 



dG 



BC 



\ 8X C dX 1 



dX L 



and 



and plugging in (4.4), we obtain 



dGAl _ qAC c bd 9Gbd 



de 

ti- 
de 



dc 

A 3 dp 

,J AB{e) - 2dXB- 



(4.11) 

(4.12) 
(4.13) 
(4.14) 
(4.15) 
(4.16) 



With these results, the linearized balance of linear momentum (4.6) becomes 

(A aA b B U b , B ) A + (B aACD G CD 0) 



3 df3 nR 

P aB = 0. 



> A 2dX B 

Assuming that A and IB are independent of X, the linearized equilibrium equations are simplified to read 



(4.17) 



J^aA^ B 



d 2 U b 



dX A dX B 

If the initial configuration is stress-free, we have 

d 2 n b 

j^aA.B ° U 



aACD G CD ^r + iJLp aB = 0. 



dX A 2dX B ' 



dX A dX B 



= -B 



aACD 



Gc 



D 



d/3 

dX A ' 



(4.18) 



(4.19) 



Let us next show that these results agree with those of classical thcrmoclasticity. We first consider a special 
class of isotropic materials. 
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Saint- Venant-KirchhofF materials. Saint- Venant-Kirchhoff materials have a constitutive relation that is 
analogous to the linear isortropic materials, namely, the second Piola-Kirchhoff stress S is given in terms of the 
Lagrangian strain E = i(C — G) as [29] 

S = A(tr E)G _1 + 2/zE, (4.20) 

or in components 

S CD = XE AB G AB G CD + 2»E CD = ^(CabG ab - 3)G CD + ^C AB G AC G BD - G CD ), (4.21) 

where A = A(X) and fi = /i(X) are two scalars characterizing the material properties. This means that S is a 
linear function of E. We next show that for this class of materials linearization of our geometric theory leads to 
linear governing equations that are identical to those of the classical linear theory of thermal stresses for linear, 
isotropic materials. We will show later in this section that this is true for any elastic material. 
We can obtain the tensor ]Q aCAB from S as follows 

,aCAB_ d f„ ab dlP \ dP aC 77tq dS CD 



B <^« Q6 _ £_ =_ — = F a D — — . (4.22) 
0G AB V dF c J 9G A b OGab 

Using 

= -G AM G BN , (4.23) 

we obtain 

B aACD G CD = -2C M NF a B (XG AB G AIN + 2 i iG AM G BN ) + (3A + 2 f i)F a B G AB . (4.24) 

The initial metric is Euclidean; in Cartesian coordinates, Gab = 5ab- Since the ambient space is also Euclidean, 
we can choose a Cartesian coordinate system whose axes coincide with the initial location of the material points 
along the material Cartesian axis. This will give, F a a = 8^, where a and A both range over 1, 2, 3. This gives 

B aACD G CD = -^±^1 5 aA . (4.25) 
Similarly, for an initially stress-free material manifold, we obtain 

A aA b B = F a MF c Ngbc \^qAM qBN + ^qABqMN + qANqBM^ (4 2g) 

For the case of an initially Euclidean material manifold with Cartesian coordinates we have 

AaA " B dxZxB = (A + fi)Ub - ab + ^ Ua ' bb - (4 ' 27) 

Therefore, Eq. (4.19) reads 

(A + /i)t/ M6 + /iC/ 0)i)6 = ^-^l^- ( 4 - 28 ) 

2 dx a 

Recalling J3 = 2a^| e , and assuming for simplicity that a is constant, we have 



dp n d dT 
2a 



dx a dx a de 

Hence 



=o 



(4.29) 



d dT 

(A + n)U bMb + fiU aM = (3A + 2fi)a— — . (4.30) 

oxa de 



=0 



In the classical theory of thermal stresses, stress-strain relations in the presence of temperatures changes can 
be written as 

(Tij = C ijk i(e k i - aS k iAT), (4-31) 
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where Cijki is the elasticity tensor. In this sense thermal strains are understood as "eigen strains" . Equilibrium 
equations in the absence of body forces read 

DAT 

^ijkl^kLj = QjfcfcQ ^— , (4.32) 

where once again we have assumed that the elasticity tensor and the coefficient of thermal expansion are 
constants. When the material is isotropic, Cijki — ^{SikSji + SuSjk) + \5ij5u and hence 

C, /./(/,/,. = //</,.,, + (A + fJt)Uj,ji, (4.33) 

and 



9 AT ,„ 9 AT 

Thus, equilibrium equations read 



^~ijkk& o 

- = (2^ + 3A)a-^— . (4.34) 



fJ>Ui t jj + (A + n)uj t ji = (2fj, + 3A)a— — . (4.35) 

Recalling that ^-| £ _ is the linearized version of temperature change, i.e., ST as e^-| e _ and that U is the 
linearized version of displacement, it is seen that the linearization of the geometric theory for Saint- Venant- 
Kirchhoff materials results in governing equations that are identical to those of the classical isotropic linear 
theory 

Let us now see if this holds for more general constitutive equations of the geometric theory. For a stress-free 
(Euclidean) initial configuration, the linearized balance of linear momentum is given by 

(A aA 6 B U b B + B aACD 8 C D 0) A = 0. (4.36) 

Assuming Gab = Sab and F a a = 5 A as above, an identity proven in [29] becomes 

A aA b B = 2<C CADB F c D F a c g cb = 2€ CADB 8 c D 8 a c 5 bc = 2€ aAbB . (4.37) 

Noting that j3 = 2a^| e=0 , (4.36) becomes 



<C aAbB U b B + B aAcc a ^~ 



flT 
1~( 



e=0 



= 0. (4.38) 



Identifying superscripts and subscripts, identifying the spatial and material indices (by aligning the material 
and spatial Cartesian coordinates as above), and using symmetries of (D, we can write 



This is identical to (4.31) if 

IBzjfcfc = — <Eijki8ki = —^ijkk- (4.40) 

Let us first show that this relation always holds for isotropic materials. For isotropic materials it can be 
shown that [25; 26; 42] 

C + ?77-G = 0. (4.41) 



Or in components 



Note that 



DC dG 

^t^Ccb + ^t^Gcb = 0. (4.42) 
dC A c oGac 

m dP AS nT1 d 2 ^ 

= <9G = <9G = (4 ' 43) 
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Differentiating (4.41) with respect to C and noting that the initial configuration is stress free we obtain 

d 2 ^ „ a 2 * 



dCd C- C+ dCdG- G = °- ^ 
Or 

ic-C + iF" 1 IB-G = 0. (4.45) 

Noting that F a A = S A and Cab = Gab, (4.45) is identical to (4.40). 

We now show that Eq.(4.40) holds even for anisotropic elastic solids. For showing this we use the fact 
that if the material is homogeneous (i.e., if the thermal expanstion properties do not depend on position) a 
uniform temperature change AT does not lead to any thermal stresses. Starting from a stress-free Euclidean 
configuration, for a uniform temperature change, one has 

SF a A = V a ,A = aAT5 a A and SS = 0. (4.46) 

We also know S = S(C, G), thus 

On on 

SS = — ■ 5C + — ■ SG = <D • SC + 2aAT F _1 B G = 0. (4.47) 

oG uG 

But note that [43] 

SCab = 9abF a A U b \B + 9abF b B U a \A = S ab S A (aAT5 b B ) + S ab S b B (aAT8 a A ) = 2aAT5 AB . (4.48) 
Substituting (4.48) into (4.47), one obtains (4.40)! In summary, we have proved the following proposition. 

Proposition. Linearization of the present geometric theory yields governing equations that are identical to 
those of classical linear elasticity. 



5 Conclusions 

In this paper, we presented a geometric theory of thermal stresses in which the material manifold is temperature 
dependent. Given a temperature distribution, the material metric is a Riemannian metric that is obtained by 
a (non-uniform) rescaling of a reference metric. In particular, starting from a Euclidean stress-free reference 
manifold, a non-uniform temperature distribution leads to a non-Euclidean material manifold. We studied the 
stress-free temperature distributions by looking at conditions that guarantee flatness of a Riemannian metric. 
We recovered some known facts from the linear theory of thermal stresses and obtained some new results for 
finite deformations. We showed that, in addition to uniform temperature distributions, there are other zero- 
stress temperature distributions. We obtained all such temperature distributions. We also studied the inverse 
problem, i.e., given a temperature distribution, what inhomogeneous coefficients of thermal expansion give zero 
stresses. In the present theory, there is no need to introduce an "intermediate" configuration. We made an 
explicit connection between our geometric theory and the previous works on multiplicative decomposition of 
deformation gradient in the presence of temperature changes. Given a temperature distribution, we obtained 
the temperature-dependent governing equations. In order to demonstrate the power of the geometric theory, we 
solved the example of an axisymmetric temperature distribution and obtained some exact results. We showed 
that linearization of the present geometric theory about a stress- free configuration results in governing equations 
that are identical to those of the classical linear thermoelasticity. 

Geometric formulation of the coupled problem of elastic deformations with heat conduction will be studied 
in a future communication. The ideas presented in this paper can also be used in modeling bodies with growing 
mass. Growth and remodeling in biological systems is an important phenomenon and a geometric study will 
shed light on the coupling between growth/remodeling and elastic deformations. 
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A Differential Geometry and Classical Geometric Elasticity 

In this section, in order to make the paper self-contained, we review some notation from geometric elasticity. For 
more details refer to [29; 1; 30]. By classical geometric elasticity we mean elasticity of bodies with stationary 
defects (if any) and a fixed material manifold. We extended this theory for thermal deformations in §3. 

For a smooth n-manifold M, the tangent space to M at a point p G M is denoted T P M and the whole 
tangent bundle is denoted TM. We denote by B a reference manifold for our body and by S the space in which 
the body moves. We assume that B and S are Riemannian manifolds with metrics G and g, respectively. Local 
coordinates on B are denoted by {A A } and those on S by {x a }. 

A deformation of the body is a C 1 embedding p : B — > S. The tangent map of p is denoted F = Tip : 
TB — > TS, which is often called the deformation gradient. In local charts on B and S, the tangent map of p is 
given by the Jacobian matrix of partial derivatives of the components of tp, as 

F = Tip : TB — > TS, Ty(X, Y) = (p(X),Dp(X) • Y). (A.l) 

If Y is a vector field on B, then p*Y = Tp • Y o ip , or using the F notation. p*Y = F • Y o p^ 1 is a 
vector field on p(B) called the push-forward of Y by p. Similarly, if y is a vector field on p{B) C S, then 
p*y = T(p^ 1 ) ■ y o p is a vector field on £> and is called the pull-back of y by p. 

The cotangent bundle of a manifold M is denoted T*M and the fiber at a point p € M (the vector space of 
one-forms at p) is denoted by T*M. If /3 is a one-form on 5, i.e., a section of the cotangent bundle T*S, then 
the one-form on B defined as 

V*/3)x • V x = /Vx) • (TV ' V X ) = /3 v( x) • (F • V x ) (A.2) 

for X G 6 and Vx G Tx£>, is called the pull-back of /3 by Similarly, the push-forward of a one-form a on 
2? is the one form on p(B) defined by p*a = (p^ 1 )*a. 

We can associate a vector field /?" to a one-form /3 on a Riemannian manifold M through the equation 

</4c,Vx> = (4,Vx) x , (A.3) 

where ( , ) denotes the natural pairing between the one form /3 X G T*M and the vector v x G T X AT and where 
(/3 X , v x} x denotes the inner product between /3 X G T X M and v x G T X M induced by the metric g. In coordinates, 
the components of 0* are given by (3 a = g ab Pb- 

(m \ 
J -tensor at X G B is a multilinear map 

T : T^B x ... x T x 2? x T x £ x ... x T x £ -> M. (A.4) 

m copies n copies 

T is said to be contravariant of order m and covariant of order n. In a local coordinate chart 

tu, : n'\v, v (! i r ...,.o] : ..w;> \ / . . . \ ;/ . (a.s) 

where a fc G T^-S and V fc G T X B. 

A two-point tensor T of type I I atXeB over a map p : B — > 5 is a multilinear map 

T : T X B x ... x T X S x TyJS x ... x T X S x T*S x ... x T*S x T x 5 x ... x T x 5 -> R, (A.6) 

•■ v ' V v ' >• v ' 11 v ' 

m copies n copies r copies s copies 

where x = <^(X). 

Let y be a vector field on S and y : B — > 5 a regular and orientation preserving C 1 map. The Piola 
transform of y is defined as 

Y = Jp*y, (A.7) 
where J is the Jacobian of p. If Y is the Piola transform of y, then the Piola identity holds: 

Div Y = J(div y) o tp. (A.8) 
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A p-form on a manifold M is a skew-symmetric [J -tensor. The space of p-forms on M is denoted by 



£l p (M). If ip : M — * N is a regular and orientation preserving C 1 map and a £ QP(f{M)) : then 

a= I f*a. (A.9) 

Let 7r : £7 — > 5 be a vector bundle over a manifold 5 and let £(<S) be the space of smooth sections of E 
and X(S) the space of vector fields on S. A connection on £7 is a map V : x E(S) — > £"(5) such that 

V f,fi,f 2 eC°°(S), Vai,a 2 eK 

*') V /lXl+/2 x 2 Y = AV Xl Y + / 2 V X2 Y, (A.10) 
ii) Vx(aiY 1 + a2 Y 2 ) = a 1 V x (Y 1 ) + a 2 Vx(Y 2 ), (A.ll) 
Hi) V X (/Y)=/V X Y + (X/)Y. (A.12) 

A linear connection on 5 is a connection on TS, i.e., V : <-f (5) x X(S) — ► X{S). In a local chart 

VaA- = 7&0*, (A.13) 

where 7^ are Christoffel symbols of the connection and di = -J~. A linear connection is said to be compatible 
with the metric of the manifold if 

V x ((Y, Z)) = {(V X Y, Z» + «Y, V X Z}) . (A. 14) 

It can be shown that V is compatible with g if and only if Vg = 0. Torsion of a connection is defined as 

T(X,Y) = VxY-V Y X-[X,Y], (A.15) 

where 

[X,Y](F) =X(Y(F)) -Y(X(F)) VFeC*°°(S), (A.16) 
is the commutator of X and Y. V is symmetric if it is torsion-free, i.e. 

V X Y- V Y X= [X,Y]. (A.17) 

It can be shown that on any Riemannian manifold (5, g) there is a unique linear connection V that is compatible 
with g and is torsion-free with the following Christoffel symbols 

^~2 9 {-d^ + d^-^J- (A - 18) 

This is the Fundamental Lemma of Riemannian Geometry [21] and this connection is called the Levi- 
Civita connection. 

Curvature tensor JZ of a Riemannian manifold (S, g) is a f ^ -tensor JZ : T*5 x T X S x T X S x T X S — > R 
defined as 

7l(a,wi,W2,w 3 ) = a (V Wl V W2 w 3 - V W2 V wl w 3 - V[ WljW2 ]W 3 ) (A. 19) 

for a €E T£S, Wi,w 2 ,w"3 G T X S. In a coordinate chart {x a } 

_,a _,e „,a _,e / a o^^ 

/<• &cd = -7T^ ~ ~Q^d + 7ce7&d - 7de7fcc- (A.20) 

Note that for an arbitrary vector field w 

™V - w a lcb = K a bcd w d + T d cb w a ld . (A.21) 

An n-dimensional Riemannian manifold is flat if it is isometric to Euclidean space. A Riemannian manifold is 
flat if and only if its curvature tensor vanishes. A Riemannian manifold (£>, G) is conformally flat if there exists 
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a smooth map / : B — > M such that G = fS, where 8 is the Euclidean metric. In isothermal coordinates the 
Riemannian metric has the following local form 

G = f(X)(dXl + ... + dXl) . (A.22) 

It is know that [2] any two-dimensional Riemannian manifold is conformally flat and the map / is unique. A 
corollary of this theorem in our theory of thermal stresses is the following. Given any smooth curved 2D solid 
that is stress free, there exists a unique change of temperature distribution such that in the new temperature 
distribution the 2D solid is flat and still stress free. Equivalently, starting from a stress free flat sheet, it is 
always possible to deform it to any smooth curved shape by changing temperature without imposing any residual 
stresses. 

For a Riemannian manifold (£>, G) the Weyl-Schouten tensor is defined as [35] 

C IJK = VkIIu - VjK IK - - A [gu^k - G ''<§Yj) ' (A - 23) 

where 7L and 1Z are the Ricci tensor and the scalar curvature, respectively. A necessary and sufficient condition 
for a Riemannian manifold (£>, G) to be conformally flat is C = when dim£> = 3. 



A.l Absolute Parallelizable (AP) Geometry 

In many physical problems in which deformation is coupled with other phenomena, e.g. plasticity, growth/remodeling, 
thermal expansion/contraction, etc. all one can hope to do is to locally decouple the elastic deformations from 
the inelastic deformations. This has led to many works that start from a decomposition of deformation gradient 
F = F e Fj, where F e is the elastic deformation gradient and F, is the remaining local deformation or inelastic 
deformation gradient. 

Given an inelastic deformation gradient, here a thermal deformation gradient, a vector in the tangent space 
of X £ B, i.e. W £ TxS is mapped to another vector W = FyW. Traditionally these vectors are assumed to 
lie in the tangent bundle of an "intermediate configuration." In the literature, intermediate configuration is not 
clearly defined and at first glance it seems to be more or less mysterious as was explained in §2.3. These are 
closely related to parallelizable manifolds (or absolutely parallelizable (AP) manifolds) going back to the works 
by Eisenhart [8, 9]. See also Youssef and Sid- Ahmed [44] and Wanas [41]. In an ?i-dimensional AP-manifold 
M, one starts with a field of n linearly independent vectors {E^)} that span the tangent vector at each point. 
Let us denote the components of ~E(A) by Ef^. The dual vectors, i.e. the corresponding basis vectors for the 

cotangent space are denoted by {E^* 1 } with components |. Note that 

Ei A) E( B) -^ and E^Efo = (A.24) 
One is then interested in equipping M with a connection Tj K such that the basis vectors {E^} are covariantly 

E(a)|j = 0. (A.25) 
E(a)|jk ~ ^(a)\kj = TZ 1 ljk^(a) + T L KJ^\a)\L- (A. 26) 



constant, i.e. 13 



Note that 
Thus (A.25) implies that 



n T L JK = (A.27) 
i.e., M is flat with respect to the connection Tj K . Note that 

9E_r 

E U)\ J = q^j +r J/<- E (A)- (A.28) 



Thus 



13 Equivalently, the tangent bundle is a trivial bundle, so that the associated principal bundle of linear frames has a section on 
M. 
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Hence 

v i _ F (A) dE U) _ { A) dE { j A) 

rjK-~Ej -qxk-Vj ~qxk~- (A.30) 

This is similar in form to the connections used by many authors, e.g. by Bilby, et al. [4] and Kondo [15] for 
dislocations, by Epstein and Elzanowski [10] for material inhomogeneities, and by Stojanovic, et al. [37] for 
thermal stresses. 

Looking at local charts and {U 1 } for the reference and intermediate configurations, we have 

dU 1 = (Ft) 1 A dX A . (A.31) 
(Ft) 1 a can be identified with Ef^, and hence 

r5^(^) V ( g^T J > (A-32) 

is the connection used in [37; 38]. Note that this connection is curvature free but has non- vanishing torsion. 



A. 2 Geometric Elasticity 

Let us next review a few of the basic notions of geometric continuum mechanics. A body B is identified with 
a Ricmannian manifold B and a configuration of B is a mapping (p : B — > S, where S is another Ricmannian 
manifold. The set of all configurations of B is denoted C. A motion is a curve c : R — > C; t i— > ipt in C. It is 
assumed that the body is stress free in the material manifold. 

For a fixed t, <^t(X) = <p(X, t) and for a fixed X, i^x(t) = ¥>(X, t), where X is position of material points in 
the undeformed configuration B. The material velocity is the map V t : B — ► R 3 given by 

V t (X)=V(X,t) = ^^ = ^ x (t). (A.33) 
Similarly, the material acceleration is defined by 

A t (X) = A(X,t) = = |v x (t). (A.34) 

In components 

^ = — + 7h Q c ^V c , (A.35) 

where 7^ is the Christoffcl symbol of the local coordinate chart {x a }. Note that A does not depend on the 
connection coefficients of the material manifold. 

Here it is assumed that ip t is invertible and regular. The spatial velocity of a regular motion ip t is defined 

as 

v t :^(£)^K 3 , vt = V t o<pt 1 , (A.36) 
and the spatial acceleration a t is defined as 

<9v 

a = v = — + V v v. (A.37) 
at 

In components 

aQ = ^ + ^" + 7&>V - (A - 38) 

Let (p : B — > S be a C 1 configuration of £> in 5, where B and 5 are manifolds. Recall that the deformation 
gradient is the tangent map of <p and is denoted by F = Tip. Thus, at each point X G B, it is a linear map 

F(X) : T x £ -> T V(X )5. (A.39) 
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If {x a } and {X A } are local coordinate charts on S and B, respectively, the components of F are 

F" A (X) = |^(X). (A.40) 

The deformation gradient may be viewed as a two-point tensor 

F(X) : T*S x T X B -> R; (a, V) ^ (a,T x¥ > ■ V). (A.41) 

Suppose B and 5 arc Riemannian manifolds with inner products ((, )) x and ((, )) x based at X G B and x G <S, 
respectively. Recall that the transpose of F is defined by 

F T : T x 5 - T x £, «FV, v)) x = «V, F T v)) x , (A.42) 

for all V G TxS, v G T X S. In components 

(i^ T (X)) A a = .g ah (x)F h B (X)G AB (X), (A.43) 

where g and G are metric tensors on S and B, respectively. On the other hand, the dual of F, a metric 
independent notion, is defined by 

F* (x) : T*S T X B; <F* (x) • a, W) = (a, F(X) W) , (A.44) 

for all a G X^S, W G Tx6. Considering bases e a and for S and 23, respectively, one can define the 
corresponding dual bases e a and E A . The matrix representation of F* with respect to the dual bases is the 
transpose of F a a- F and F* have the following local representations 

F = F a A — ® dX A , F* = F a A dX A ® — . (A.45) 
ox a ox a 

The right Cauchy- Green deformation tensor is defined by 

C(A) : T X B -► T x £, C(X) = F(X) T F(X). (A.46) 

In components 

C A B = (F T ) A a F a B . (A.47) 

It is straightforward to show that 

C b = p*(g) = F*gF, i.e. C AB = (g ab ° (p)F a AF b b- (A.48) 

Let <pt : £> — > 5 be a regular motion of £> in 5 and P C S a p-dimcnsional submanifold. The Transport 
Theorem says that for any p-form a on 5 

~dt 



a = L v a, (A.49) 



where v is the spatial velocity of the motion. In a special case when a = fdv and V — U is an open set, one 
can write 



d 
dt 



fdv = 



% + div(/v) 



dv. (A.50) 



